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We have investigated the center-of-mass oscillations of 
a 87 Rb Bose-Einstein condensate in an elongated magneto- 
static trap. We start from a trapped condensate and we 
transfer part of the atoms to another trapped level, by apply- 
ing a radio- frequency pulse. The new condensate is produced 
far from its equilibrium position in the magnetic potential, 
and periodically collides with the parent condensate. We dis- 
cuss how both the damping and the frequency shift of the 
oscillations are affected by the mutual interaction between 
the two condensates, in a wide range of trapping frequencies. 
The experimental data are compared with the prediction of a 
mean-field model. 



is damped and its frequency is shifted upwards, with respect 
to the non interacting case. 

The Paper is organized as follows. In Section [n|, we describe 
the experimental setup. Then, in Section [lll| , we present the 
experimental results, analyzing the effect of mutual interac- 
tions on the center-of-mass motion over a wide range of trap 
frequencies. The data are compared with the predictions of 
the Gross-Pitaevskii theory for two coupled condensates at 
zero temperature, according to the model discussed in Ref. 
|l2j . Finally, in Section ^ we draw the conclusion. 

II. EXPERIMENTAL SETUP 



I. INTRODUCTION 

The issue of the interaction between Bose-Einstein conden- 
sates in different internal states has deserved much attention 
in the recent literature, both from the experimental and 
theoretical point of view The study of these mixtures 

of quantum fluids includes several important topics, as the 
characterization of the dynamical behaviour, the response in 
presence of an external coupling, the phase coherence and the 
superfluid properties of the system. 

Recently we have demonstrated an experimental method 
for a sensitive investigation of the interaction between two 
condensates [^, |IoUll| , where two condensates are made to col- 
lide after periods of spatially separated evolution and quan- 
titative information is extracted from the resulting collective 
dynamics. 

In this Paper we report on new measures performed on 
such system in a wide range of trapping frequencies, which 
allow a systematic investigation of the effects of the collisions 
between the two condensates. 

The experimental setup follows the scheme reported in Ref. 
p. We use a radio-frequency (rf) pulse to produce two 

Rb condensates in the states \F = 2, irif = 2) = |2) and 
2, 1) = |1). Due to the different magnetic moments and the 
effect of gravity, they are trapped in two potentials whose min- 
ima are displaced along the vertical y axis by a distance much 
larger that the initial size of each condensate. As a conse- 
quence the 1 1} condensate, initially created with the same den- 
sity distribution and at the same position of the |2) conden- 
sate, undergoes large center-of-mass oscillations. Moreover, 
the two condensates periodically collide, and this strongly af- 
fects the collective excitations of both the condensates |q]. In 
particular the center-of-mass oscillation of the 1 1} condensate 



We start with a single |2) condensate containing typically 
1.5 x 10 s 87 Rb atoms at a temperature below 130 nK, confined 
in a 4-coils Ioffe-Pritchard trap elongated along the horizontal 
z symmetry axis |13[| , which produces an axially symmetric 
harmonic potential. The harmonic oscillator axial frequency 
is held fixed to v Z 2 = 12.6 Hz, while the radial frequency 
v±2 has been varied in the range 113 -f- 310 Hz by chang- 
ing the minimum magnetic field Bq according to the relation 
i/±2 = 221 Hz/Bo [G] 1//2 . The density distribution is an in- 
verted parabola of the Thomas- Fermi (TF) regime jlf|, and 
the typical radial and longitudinal radii of the condensate are 
R± = 4.3 -7- 2.4 fim and R z — 39 4- 58 /zm, in the range of 
frequencies considered here. 

To populate the j 1} state we apply a rf oscillating magnetic 
field. After the rf pulse the initial |2) condensate is put into 
a coherent superposition of different Zeeman |m/} sub- levels 
of the F — 2 state, which then move apart: |2) and |1) are 
low-field seeking states and stay trapped, |0) is untrapped and 
falls freely under gravity, while | — 1} and | — 2} are high-field 
seeking states repelled from the trap. 

By fixing the duration and varying the amplitude of the rf 
field, we control the relative population in different Zeeman 
sub-levels juj. In the experiments described here we use a 10 
cycles rf pulse at 1.24 MHz with an amplitude B=20 mG to 
quickly transfer ~ 13% of the atoms to the |1) state without 
populating the |0), | — 1} and | — 2} states. 

In the frequency domain, the pulse width exceeds by almost 
one order of magnitude the transition broadening due to the 
inhomogeneous magnetic field across the condensate. As a 
result, the coupling strength does not depend on position and 
1) is produced with the same spatial density profile as |2). 
Furthermore, since the characteristic time scale of the collec- 
tive dynamics is of the order of 1 ms, no evolution occurs for 
the wavefunctions during the pulse, but, immediately after, 
excitations start in both the condensates. 

The |1) condensate experiences a trapping potential with 
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lower axial and radial frequencies (vi = U2/V2), whose min- 
imum is displaced along the vertical y axis by a distance 
2/o = g/^±2 due to gravity ("gravitational sagging"). In the 
frequency range considered here yo varies between 19 and 
2.6 /im. After the rf-pulse, the |1) condensate moves apart 
from 1 2), and begins to oscillate around its equilibrium posi- 
tion. Due to the mutual repulsive interaction, the latter starts 
oscillating too, though with a much smaller amplitude. Actu- 
ally, each condensate moves in an effective potential which is 
the sum of the external potential (magnetic and gravitational) 
and the mean-field one. Here we have studied the dynamics 
of the |1) condensate by varying the permanence time in the 
trap up to 40 ms. 

After the release from the magnetic trap, we let the clouds 
drop and expand for 30 ms, and finally we take an absorption 
image on a CCD array. We notice that the process of switch- 
ing off the magnetic field is not instantaneous, and lasts about 
1 ms. Moreover it produces a magnetic gradient which affects 
the initial velocity of the two condensates. The acquired ve- 
locities experimentally observed are vi v = 0.7± 0.1 cm/s and 
t>2 y = 1.4 ±0.1 cm/s, respectively for the |1) and |2) conden- 
sates. In the following we use an unknown, effective, switch- 
off delay time as the only fitting parameter for our theoretical 
model. 




FIG. 1. Images of the two condensates before and after the 
collision (from left to right). 



III. CENTER-OF-MASS OSCILLATIONS 

Let us recall the basic features of our system, as reported for 
the original experiment in Ref. |(|. The initial configuration 
corresponds to the stationary ground-state with all the N 
trapped atoms in the \2) condensate. Afterwards, at t = 0, 
N\ ~ 0.13iV atoms are transferred from the \2) to the |1) 
state, the former remaining with N2 = N — Ni atoms. 

The dynamics of the condensate |2) is characterized by 
small amplitude center-of-mass oscillations and by shape os- 
cillations which are both affected by the mutual repulsion 
with the other condensate when they periodically collide (see 
Fig- 0) I- 

Here we are mainly interested in the motion of the |1) con- 
densate, which is on the contrary dominated by large oscil- 
lations, whose frequency v\ turns out to be larger than the 
"bare" trapping frequency v\_\. In fact the condensate |1) 
moves in an anharmonic effective potential (see Fig. ^ due to 
the presence of the condensate |2), which provides an extra 
repulsion with a consequent up-shift of frequency. This shift 
is one of the clear signatures of the mutual interaction, since 
it only occurs if the two condensates overlap periodically. 




Vertical position 
FIG. 2. Trapping potentials for the two condensates in the 
vertical direction for i^j_2 = 116 Hz (lengths are in units of 
a±2 = [ft/(™<-'±2)] 1/ ' 2 — 1-00 /jm). The |1) condensate moves 
in an effective potential (solid line) which is modified by the 
mean field repulsion of |2) (right side). 

Furthermore, the oscillations appear to be damped. At the 
level of the Gross-Pitaevskii theory this damping is not due 
to dissipative processes (the total energy is conserved), but 
to a redistribution of the kinetic energy initially associated 
with the center-of-mass motion of the |1) condensate, which 
is eventually shared among the other degrees of freedom of the 
system (center-of-mass of the other condensate and internal 
excited states of both of them) . 

To investigate the role of the interactions in the large 
center-of-mass oscillations of the |1) condensate in different 
dynamical regimes, we have carried out a systematic analysis 
by varying the trap frequency v± 1 in the range 80 -f- 220 Hz 
(which corresponds to the above mentioned range 113 -f- 310 
Hz for u±a). 

In Fig. ^ we show the center-of-mass motion of the |1) 
condensate, for some values of the radial trapping frequency, 
by comparing the experimental data with the prediction of 
the 2D Gross-Pitaevskii model discussed in Refs. [^,^|. In 
this case the only fitting parameter in the theoretical curves 
is the delay time caused by the non instantaneous switch-off of 
the trapping potential. Its value is ~ 0.7 ms. We notice also 
that in our model the effect of the expansion on the center-of- 
mass motion is considered as a simple free fall under gravity, 
neglecting any residual interaction between the condensates. 
Actually, although the two clouds acquire different velocities 
during the release from the trap, and cross each other in the 
fall, we expect them to quickly become so dilute, due to the 
fast radial expansion, that we can neglect their mutual mean- 
field repulsion. This is also supported by numerical simulation 
including interaction in the early stages of expansion Jl6[. 
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FIG. 4. Relative shift of the |1) oscillation frequency versus 
its "bare" trapping frequency, i.e. in absence of interactions 
with 1 2). The experimental data are compared with the pre- 
diction of the classical model (solid line) in Ref. [^] and of 
the 2D-GP model (empty triangles). 

The values for the frequency shift and the damping time in 
the whole range of frequency are summarized in Figs. ^ and 
|B] respectively. In Fig. ^ the relative shift of the 1 1) oscillation 
frequency is plotted versus its "bare" trapping frequency, i.e. 
in absence of interactions with |2). The experimental data 
are compared with the prediction of the classical model (solid 
line) and of the 2D-GP model (triangles). In Fig. ^, we show 
the "quality factor" Q — u±i ■ r of the oscillator as a function 
of t being the damping time, comparing experimental 
points and theoretical predictions for our 2D model. 
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FIG. 3. Center-of-mass oscillations of the |1) condensates 
after the ballistic expansion for three values of the "bare" 
trapping frequency (u±i = 92.1, 134.6, 219.2 Hz, from top to 
bottom). The measured center-of-mass positions (arbitrary 
units) are plotted as a function of the trapped evolution time 
(in units of w^[ 2 ), and compared with the prediction of our 
2D model (solid line). 
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FIG. 5. Measured "quality factor" Q = lu±\t for the |1) 
oscillation versus its "bare" frequency compared with the pre- 
diction of the 2D-GP model (empty triangles). 

We notice that, concerning the frequency shift, there is a 
general agreement (within the experimental uncertainty) be- 
tween GP theory and experiment in the whole range of fre- 
quency, except for the frequency of the original experiment 
(116.3 Hz). By the way, we have found that in this frequency 



range around 116.3 Hz there is an almost degeneracy between 
the energies of the two condensates (mean field + kinetic + 
potential energy that the two condensates can exchange dur- 
ing the collisions, see Ref. Jl^]). At this stage, we cannot 
tell whether this can play some physical role, or it is a pure 
coincidence. 

On the contrary the observed damping is generally under- 
estimated by the 2D GP model. This is not surprising, since 
there could be other sources of damping not included in our 
mean-field description. To be more specific, one should con- 
sider the following points. 

(i) The theoretical curves are based on a 2D model which 
only includes the radial degrees of freedom in the a;y-plane. 
In fact, being the geometry of our condensates strongly elon- 
gated in the longitudinal direction, we have approximated the 
system as a uniform cylinder, by freezing out the axial dynam- 
ics (see Ref. |l| for a comprehensive discussion). Therefore, 
the remaining contribution to the damping could in part arise 
from the excitations of axial modes (possible distortions of the 
shape of the two condensates along z), which are ignored in 
the present analysis. Anyway, there are indications from pre- 
liminary 3D simulations that these contributions should not 
be so important, at least for the case at 116.3 Hz. 

(ii) We observe a certain fragmentation of the |1) conden- 
sate after some collision. This may have an influence on the 
measured center-of-mass since the density distribution devel- 
ops low density asymmetric wings, which we might be unable 
to detect properly. 

(Hi) From the analysis of the experimental data there are 
indications in support of the occurrence of atom (and energy) 
losses from both condensates. In fact, the relative number 
of atoms for the |1) state (Ni/N) is in average a decreasing 
function of the trapping time. We remind that, although 
each experimental point corresponds to a different experiment 
(each point corresponds to a condensates released from the 
trap after a variable time), they are all performed with almost 
the "same" initial conditions. In principle we can exclude 
that these losses originate from the condensate finite lifetime 
(0.7(1) s) or heating (dT/dt = 0.11(2) (jK/s), due to the short 
permanence times in the trap (less than 40 ms). Actually, 
this behaviour could be explained with a loss of atoms due 
to binary collisions between atoms of the two condensates, an 
effect that is not included in our GP mean-field theory. 

Concerning this last point, an appealing mechanism could 
be the occurrence of elastic scattering processes between 
atoms of the two overlapping condensates, as discussed in 
Ref. jl7| for the impurities motion inside a Bose condensate 
[ p^[ and in the analysis of matter 4- waves mixing jljj. In fact, 
since Bose-Einstein condensates are superfluids, when the rel- 
ative velocity of the two condensates exceeds a critical value of 
the order of the sound velocity, we could expect the GP mean- 
field approximation to break down and "elastic scattering" 
losses to take place. Indeed, in our experiment, the two con- 
densates collide at a velocity which is larger than the (maxi- 
mum) speed of sound at the top of j2), Co = hy/^Kn^w^jm, 
by a factor a that can be tuned by changing the radial fre- 
quency. In the frequency range considered, a goes from 4.2 
at v±2 ~ 113 Hz to 2.1 at v±2 = 310 Hz. So, we should ex- 
pect the effect of the "elastic scattering", if any, to be more 
pronounced at low radial frequency. As a consequence the 



results of the GP model should be more accurate at high fre- 
quency. This conclusion, however, is not supported by the 
experimental observations. 

A detailed theoretical analysis of our experiment in this 
respect is under way at the NIST group, in Gaithersbourg 
(USA). Preliminary results give a better agreement for the 
damping time with the inclusion of elastic scatterin g |lql , at 
least for the case of the original experiment in Ref. ]12j|. 



IV. CONCLUSIONS 

To summarize, we have analyzed the large center-of-mass 
oscillations of a Bose-Einstein condensates in a strongly elon- 
gated magnetic trap, over a wide range of frequencies. The 
condensate is produced by transferring part of the atoms from 
a parent condensate in the \F = 2,m/ = 2) state to the 
\F — 2,m/ = 1) state. The new condensate is created far 
from its equilibrium position in the magnetic potential, and 
undergoes large oscillations, periodically colliding with the 
parent condensate. We have shown that the mutual inter- 
actions between the two condensates deeply affect both the 
amplitude and the frequency of the oscillations. The experi- 
mental data are compared with the prediction of a 2D model 
based on the Gross-Pitaevskii theory for two coupled conden- 
sates at zero temperature. This model accounts for the basic 
features of the system, and the agreement with the data is 
generally satisfactory. Nevertheless, the discrepancy observed 
for the damping (and in some case also for the frequency) 
suggests that processes of atom losses (e.g. elastic scattering 
losses |17| ) could play a significant role beyond the mean-field 
theory. 
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